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We have studied the driving rate and temperature dependence of the power-law exponents that 
characterize the avalanche distribution in first-order phase transitions. Measurements of acoustic 
emission in structural transitions in Cu-Zn-Al and Cu-Al-Ni are presented. We show how the 
observed behaviour emerges within a general framework of competing time scales of avalanche 
relajcation, driving rate, and thermal fluctuations. We have confirmed our findings by numerical 
simulations of a prototype model. 
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The dynamics of first-order phase transitions (FOPT) 
in the presence of disorder is a long-standing problem 
which is undergoing renewed interest. Besides the phe- 
nomena of nuclcation, metastability and hysteresis, ad- 
ditional features arc presently being revisited: (i) the 
extension of the transition over a broad region on the 
generic temperature-field phase diagram (where the field 
is the intensive parameter conjugated to the order pa- 
rameter) ; and (ii) the fact that transitions occur through 
a sequence of avalanches that link metastablc states when 
thermal fluctuations arc not dominant at low enough 
temperatures. These order parameter discontinuities are 
associated with the sudden transformation of a fraction of 
the system. This behaviour has been observed in many 
magnetic Q, Q, S 13 , superconducting [5|, |a| and struc- 
tural |3, |8J transitions and in vapour condensation on 
porous media [9j. In many cases this phenomenon has 
been described as "athcrmal" FOPT J8|. In practice, this 
means that for the FOPT to proceed it is necessary to 
drive the system externally by varying the temperature 
T or the generic field H . When T and H are constant no 
avalanches are detected even in the case of long waiting 
times. This reflects the fact that thermal fluctuations are 
irrelevant compared with the high-energy barriers sepa- 
rating the metastablc states, and that temperature acts 
as a scalar field. The transformed fraction as a function 
of time only depends on the present and past values of T 
or H , but not explicitly on time. 

In different systems, avalanches have been detected 
with a variety of experimental techniques inclu ding 
induction (Barkhausen) |j, llOj, magnetization |lll| . 
calorimetry [ij], resistivity [5j, acoustic emission (AE) 
Q , capacitance [3 and optical measurements [l|, y, llJl ■ 
Interestingly, in many cases the statistical distribution of 
the avalanche properties (area, duration, released energy, 
signal amplitude, etc..) exhibits power-law behaviour 
over several decades. This has suggested the existence 
of criticality and has stimulated the search for univer- 
sality in the power-law exponents. However, many dif- 
ferent exponent values have been reported and compar- 



ison with the theoretically predicted universality classes 
[ij, list llq is controversial. Experimentally accessi- 
ble quantities such as the distribution of signal ampli- 
tudes are not straightforwardly comparable with those 
obtained in numerical simulations such as the avalanche 
sizes (which arc well defined). Moreover, the experimen- 
tal distributions in some cases show exponential correc- 
tions which make it difficult to obtain the power-law ex- 
ponents. 

Systematic measurements with well-controlled exter- 
nal parameters are scarce. In structural transitions 
the exponent associated with AE distributions has been 
shown to depend on repetitive cycling through the transi- 
tion. Stationary values, which can be grouped into differ- 
ent classes, are obtained after a relatively large number 
of cycles [l3|. For Barkhausen measurements in Fe-films, 
an increase of the exponent has been found when tem- 
perature is decreased from room temperature to 10 K 
[iSj . For soft-magnetic materials two classes have been 
proposed [13: one exhibiting exponents which decrease 
with increasing driving rate r and the other with expo- 
nents which are independent of r |2f]| . 

In this Letter we present systematic measurements of 
the influence of cooling rate r = |r| on the distribu- 
tion of amplitudes of AE signals in structural transitions. 
We analyze two different systems with different degree of 
athermal character. In one case, the exponent has been 
found to decrease with r, while it increases in the other 
case. Results are well explained within a general frame- 
work, which accounts for the effect of thermal fluctua- 
tions. Numerical simulations corroborate the proposed 
framework. 

The studied samples are Cugs.oZnig.oAlie.o and 
Cu68.4Al27.8Ni3. 8 single crystals which were investigated 
in a previous work [3. On cooling, these alloys trans- 
form martensitically at Tm = 245 K and Tm = 257 
K, respectively. Both transitions extend in temperature 
AT '^ 35 K. Previous results |31 conflrmed that the Cu- 
Zn-Al behaves athermally for temperature driving rates 
larger than a value which is less than 0.1 K/min. This 



was verified by checking that the AE activity per tem- 
perature interval as a function of T is independent of r 
(scahng). In contrast, Cu-Al-Ni is much less "athermal". 
Although it clearly shows avalanche dynamics, the activ- 
ity does not scale at the studied driving rates and few 
isothermal signals occur. It was estimated that scaling 
for this sample will hold for r above 50 K/min. 

The AE signals were detected by a piezo-electric trans- 
ducer using the experimental set-up explained elsewhere 
[SJ . The samples were cycled through the transition more 
than 100 times in order to reach the stationary path. In- 
dividual AE signals were detected during cooling runs 
at selected rates and their amplitudes A were analyzed. 
The experimentally accessible range is 4 x 10~^V< A < 
1.5 X 10~'^V. For the Cu-Al-Ni sample averages over ^ 10 
cycles were also performed in order to obtain good statis- 
tics. The typical number of analyzed signals for each 
rate is ~ 10^. The temperatures between which the cy- 
cles are performed are always the same: T = 220 K and 
T = 340 K for the Cu-Zn-Al sample and T = 220 K and 
T = 320 K for Cu-Al-Ni. The distribution of AE ampli- 
tudes (as well as the size distribution in the simulations 
presented below) is studied by the maximum likelihood 
method which is independent of the way histograms are 
plotted. We have fitted the data to the normalized prob- 
ability law: p{A) ~ A~'^e~'^^. Two different fits are 
performed: first, a one-parameter fit of the exponent a 
by imposing A = and, second, a two-parameter fit of 
a and A. The values of a reported here correspond to 
the first fitted value, but error bars include the second 
fitted value. For the second fit |A| < 1000 V'^ in aU 
cases, which indicates that differences between the two 
fitted functions start to be important outside of the ex- 
perimental window. Figure ^ shows the dependence of a 
on the rate r. The striking feature is that for the most 
athermal case (Cu-Zn-Al) a decreases with r while for 
the less athermal sample (Cu-Al-Ni), it increases. More- 
over, in both cases, the dependence on r is logarithmic, 
as shown by the fitted lines. A decrease of the exponent 
of Barkhauscn avalanche sizes (area below the signals) in 
polycrystalline Si-Fe with increasing rate is reported in 
Ref. |2fl|. We have plotted these data on a log- linear scale 
in the inset of Fig. ^ and a good logarithmic dependence 
is also observed. 

Understanding of the results comes from the analy- 
sis of the relevant time scales involved in the problem. 
These are: (i) The time scale of the avalanche relaxation 
Tav ■ Although avalanche durations are known to be also 
power-law distributed over several decades, the detected 
values range between Tav ~ 10~^ — 10"'^ s. (ii) The char- 
acteristic time associated with the driving rate Tdr ■ This 
can be conveniently defined as a function of the extension 
of the FOPT as r^r = AT/r. In our case from the ex- 
perimentally accessible rates Tdr ^ 10^ — 10^ s. (iii) The 
third time scale is associated with the activation of the 
FOPT by thermal fluctuations Tth- This scale decreases 
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FIG. 1: Fitted exponents corresponding to the amplitude dis- 
tribution of AE signals as a function of the driving rate r . The 
inset shows the exponents corresponding to the Barkhausen 
signal integrated amplitudes from Ref. |23 



with temperature and depends on the free energy barri- 
ers. For our samples, given the above-mentioned limiting 
rates for the athermal behaviour, we have estimated that 
Tth > 10^ s for Cu-Zn-Al and Tt,, < 10^ s for Cu-Al-Ni. 
Moreover, given that AT/Tm ^ 0.1 is small for the stud- 
ied samples, we expect that, in both cases, Tth does not 
depend on r. 

The dynamics of the FOPT will therefore be deter- 
mined by competition between these three time scales. In 
the "adiabatic" limit Tdr/Tth ~* and Tav/Tdr ~* 0, i.e., 
when the time scales are well separated Tav ^ Tdr ^ Tt/i, 
the avalanches are well defined and both temperature 
and finite driving rate do not affect avalanche scaling. 
When Tdr approaches Tav, the system will still display 
avalanche behaviour but, will start to overlap. Typ ically, 
small avalanches merge to form larger signals |21|. This 
leads to a decrease in the exponent a when Tav/Tdr on r 
increases. On the other hand, by increasing tempera- 
ture and or decreasing driving rate, competition of the 
thermal fluctuations and driving may produce diverse ef- 
fects. Strictly speaking, for any nonzero value of the ratio 
Tdr/Tth oc 1/r the activity may not stop even at r = 
until thermal equilibrium is reached, making the concept 
of avalanches difficult to define. Nevertheless, from the 
experimental point of view, there is always a threshold 
below which signals cannot be detected. Thus, distinct 
signals can still be measured for low enough values of 
Tdr/Tth- We expect that increasing thermal activation 
(decreasing Tth) will help the advance of the FOPT by 
promoting the transformation of domains which other- 
wise would not transform under athermal circumstances, 
thus resulting again in the merging of small avalanches. 
On the other hand, thermal fluctuations may initiate 



small avalanches when the system is slowly driven (r^r 
is large) which will not initiate at such early fields when 
the system is driven fast. The potential effects of the 
two time scale ratios TavJTdr and Tdr/Tth on the scaling 
exponents are shown schematically in Fig. 13 Generally, 
for a given system at finite temperature, driving rate de- 
pendence is determined along the projected curve on the 
basal plane. Thus, depending on the experimentally ac- 
cessible range of driving rates we can find an increase of 
the exponent for the less athermal samples, a region of 
constant exponents, or a decrease of the exponent for the 
more athermal samples. This enlarged space of parame- 
ters offers a suitable scenario to understand the present 
results for the acoustic emission measurements, as well 
as the behavior observed by large temperature changes 
recently reporte d l lSl . The studied case of Barkhausen 
avalanches |l9l 1201 represents the limit where thermal 
fluctuations can be neglected. 
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FIG. 2: Schematic diagram showing the behaviour of the 
exponent as a function of the relevant time scales. The hy- 
perbola on the base plane shows the behaviour at constant 
temperature when increasing the driving rate. 

In order to substantiate the above arguments we have 
performed numerical simulations. The prototype model 
for hysteresis and avalanches in athermal FOPT is the 
3D-Random Field Ising Model (RFIM) at T = with 
metastable dynamics [ij . It consists of a ferromagnetic 
Ising model with local random fields which are Gaussian 
distributed with standard deviation a. The system is 
driven from the saturated state by decreasing the exter- 
nal field H . Spins flip according to a local relaxation 
rule which is rcponsiblc for hysteresis. By construction, 
the model is in the perfect adiabatic limit Tav ^ Tdr ^ 
Tth — oo. The flrst inequality is guaranteed because as 
soon as one avalanche starts, H is kept constant until 
the system reaches a new metastable state. We will refer 
to such dynamics as athermal-adiabatic dynamics. The 
model reproduces the existence of avalanches whose size 
s and duration At can be statistically analyzed. In a 
large region close to the critical point ac — 2.21 |23 |. 
the avalanche sizes corresponding to a full half-loop dis- 
tribute, to a good approximation, according to a power- 
law. Sensitivity of the exponent to system size L and a 
has recently been studied in detail |22, I23. Moreover, 



given that the distributions are only approximate power- 
laws, the exponents also depend on the fltting range. 

Extension of the model in order to incorporate the fl- 
nite driving rate r = \H\ and thermal fluctuations ef- 
fects is not straightforward. The appropriate enlarge- 
ment of the parameter space can be achieved in several 
ways. Here we adapt the approach in which the sepa- 
ration of signals is preserved to a large degree. A first 
attempt to incorporate the finite driving rate, which we 
refer to athermal- step- driven dynamics, was proposed in 
Refs. |2J, |2^ (Here we extend this work to lower driv- 
ing rates). It consists of increasing the external field 
H hy a certain step AiJ and keeping it fixed until a 
new metastable situation is reached. The obtained sig- 
nals are superpositions of the avalanches of the athermal- 
adiabatic model. The statistical analysis shows that the 
exponent decreases with the increasing "rate" (quantified 
by /S.H). Results for a system with a = 2.0 are shown in 
Fig. O as filled circles |28| . Data are compared with the 
exponent corresponding to athermal-adiabatic dynamics 
which is shown by a dashed horizontal line. A logarithmic 
decrease with r is obtained in a broad region, in quali- 
tative agreement with the results obtained for Cu-Zn-Al 
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FIG. 3: Exponents corresponding to the avalanche size dis- 
tribution as a function of the driving rate obtained by dif- 
ferent numerical algorithms at different temperatures T, as 
explained in the text. 

Strictly speaking, for the athermal-step-driven dynam- 
ics, r is not constant along the hysteresis path. This 
could be important since long avalanches are known to 
concentrate on the central part of the hysteresis loop and 
correlations may exist. In order to check whether or not 
this affects the results, we propose other dynamics which 
we will call athermal- finite- driving-rate dynamics. Us- 
ing the original athermal-adiabatic dynamics, one records 
the sequence of avalanches in half a loop (including the 
field value for which they occur, their size, and their du- 



ration). For a given finite rate H, it is possible to analyze 
the sequence of recorded avalanches and determine which 
of them will overlap. Thus, one can obtain the modified 
sequence of signals and perform the statistical analysis. 
Results obtained by this method are indicated in Fig. |2| 
as open circles. In order to compare it with the previ- 
ous data, H has been multiplied by the average duration 
(At) of the pulses in a half loop. Agreement with the 
previous dynamics is very good, thus supporting the use 
of step-driven dynamics in simulations. 

The numerical results obtained arc in agreement with 
the behaviour proposed in Fig. [5| and correspond to the 
Tdr/Tth = plane. To simulate the effect of a finite Tth 
we must modify the dynamics proposed above introduc- 
ing the effect of thermal fluctuations, but still keeping 
separate signals. This is achieved by slightly modifying 
the step-driven algorithm. After each AH, not only un- 
stable spins relax, but also an extra small fraction of lo- 
cally stable spins is reversed towards the new phase due 
to thermal activation. This is done "a la Metropolis", 
simultaneously and independently for each spin, with a 
probability p = min(l, e~ ''^), where AE is the energy 
change associated with the reversal of each spin. Dur- 
ing the subsequent evolution of the avalanche, the field is 
kept constant and temperature is set to 0. The exponents 
for three different temperatures are shown in Fig. 13 As 
can be seen non-monotonous behaviour is obtained, as 
qualitatively illustrated in Fig. [3 

Several models have been previously proposed for the 

understanding of the influence of driving rate on the 

ower-law exponents characterizing signal distributions 

I21L |2y, 1231 ■ All of them correspond to athermal situ- 
ations and account for the decrease of the exponent with 
increasing r, as a consequence of the overlap of small 
avalanches. In some cases a linear decrease of the ex- 
ponents with r is predicted. Our experiments and sim- 
ulations fit better with a logr dependence. Moreover, 
we have also shown that Barkhausen data in Ref. |20j 
also conform to this logarithmic tendency. We expect 
that the origin of this functional dependence lies in the 
correlation between avalanches evolving in the restricted 
geometry of the transforming sample rather than in the 
dynamics itself. Such correlations have not been taken 
into account in previous models which correspond to sin- 
gle advancing interfaces in magnets [J, |2a, |23| or systems 
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with uncorrclatcd avalanches 

In conclusion, we have shown that the avalanches at 
FOPT may exhibit a variety of the driving rate depen- 
dences related to the interplay between the three relevant 
time scales. We have, therefore, clarified under which ex- 
perimental circumstances one can expect exponents sim- 
ilar to those predicted in the athermal adiabatic limit. 
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The actual values of the exponent obtained with the dif- 
ferent numerical algorithms discussed in this work de- 
pend on a, L and the fitting range of s. We have verified 
that for 2.0 < a < 2.4 and L > 20, the qualitative be- 
haviour for 2 < s < 1000 is that shown in Fig.|H] 
The actual values of the exponents are not comparable 
since the simulations correspond to avalanche sizes and 
the experimental data correspond to amplitudes of AE 



signals. 



